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We study the collapse in a coupled Bose-Einstein condensate of two types of bosons 1 and 2 under 
the action of a trap using the time-dependent Gross-Pitaevskii equation. The system may undergo 
collapse when one, two or three of the scattering lengths aij for scattering of boson i with j, i,j — 1,2, 
are negative representing an attractive interaction. Depending on the parameters of the problem a 
single or both components of the condensate may experience collapse. 
PACS Number(s): 03.75. Fi,05.30.Jp 



I. INTRODUCTION 

The experimental detection of Bose-Einstein con- 
densation (BEC) at ultralow temperature in dilute 
trapped bosons (alkali metal and hydrogen atoms and 
the recent possibility in molecules) has spurred intense 
theoretical activities on various aspects of the conden- 
sate Many properties of the condensate are usu- 
ally described by the mean-field time-dependent Gross- 
Pitaevskii (GP) equation |^,||. One of the most inter- 
esting features of BEC has been observed in the case of 
attractive interatomic interaction In that case the 
condensate is stable for a maximum critical number of 
atoms, beyond which the condensate experiences a col- 
lapse. When the number of atoms increases beyond the 
critical number, due to interatomic attraction the radius 
tends to zero and the central density of the condensate 
tends to infinity. Consequently, the condensate collapses 
emitting particles until the number of atoms is reduced 
below the critical number and a stable configuration is 
reached. The condensate may experience a series of col- 
lapses [§|J3|. This phenomenon was observed in the BEC 
of ^Li atoms with negative scattering length denoting at- 
tractive interaction where the critical number of atoms 
was about 1400 1^,^. Theoretical analyses based on the 
GP equation in the case of ^Li atoms also confirmed this 
collapse illi. 

More recently, there has been experimental realiza- 
tion of BEC involving atoms in two different quantum 
states 1^,0. In one experiment *^Rb atoms formed in 
the F = 1, m = —1 and F = 2, m — 1 states by the 
use of a laser served as two different species, where F 
and m are the total angular momentum and its pro- 
jection 1^. In another experiment a coupled BEC was 
formed with the ^^Rb atoms in the F = 1, m = — 1 and 
F = 2, m = 2 states 10|. It is possible to use the same 
magnetic trap to confine atoms in two magnetic states 
and this makes these experimental investigations techni- 
cally simpler compared to a realization of BEC with two 
different types of atoms requiring two different trapping 
mechanisms. This is why so far it has not been possi- 
ble to prepare a coupled BEC with two different types 
of atoms. In addition to coupled atomic condensates. 



there has been consideration of a hybrid BEC where one 
type of bosons are atoms and the other molecules [pd]| . 
These initiated theoretical activities in BEC involving 
more than one types of bosons using the coupled GP 
equation 11 1^. 

In addition to just forming a coupled BEC with two 
quantum states of the same atom, these studies also 
yielded crucial information about the interaction among 
component atoms and measured the percentage of each 
quantum states in the condensate [^,|0|. It has been 
found that ^''Rb atoms have repulsive interaction in all 
three quantum states. Also, the strength of repulsive in- 
teractions inF=l,m = — 1 and F — 2,m — 2 states are 
essentially identical. The interaction between an atom in 
the F — 1, m = — 1 state and another in the F = 2, to = 2 
state is repulsive. As the change in the m value of a 
atomic quantum state does not correspond to a substan- 
tial structural change, it is likely that such change would 
not correspond to a large change in the atomic interac- 
tion. 

Here we study theoretically the collapse in a cou- 
pled BEC composed of two quantum states 1 and 2 of 
a bosonic atom using the coupled time-dependent GP 
equation. We motivate this study by considering two 
possible atomic states of ^Li whenever possible. An ex- 
periment of collapse in a coupled BEC has not yet been 
realized but could be possible in the future. In the case 
of ^Li the interaction in state 1 is taken to be attractive 
which is responsible for collapse. Here there are three 
types of interactions denoted by the scattering lengths 
o-yi ijj = lj2, between states i and j. A negative (pos- 
itive) scattering length denotes an attractive (repulsive) 
interaction. We study the collapse with different possibil- 
ities of attraction and repulsion between atoms in state 
1 and 2. If one of the scattering lengths is negative, at 
least one component of the condensate may experience 
collapse. If two of the scattering lengths are negative 
one can have collapse in both components. Specifically, 
one can also have collapse of both components if ai2 is 
negative and an, i = 1,2 are positive. 

The usual GP equation conserves the number of atoms. 
The dynamics of the collapse (growth and decay of num- 
ber of atoms) is best studied by introducing an absorp- 
tive contact interaction in the GP equation which allows 
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for a growth in the particle number from an external 
source. One has also to introduce an imaginary quartic 
three-body interaction term responsible for recombina- 
tion loss from the condensate |||. If the strengths of 
these two terms are properly chosen, the solution of the 
time-dependent GP equation could produce a growth of 
the condensate with time when the number of atoms is 
less than the critical number. Once it increases past the 
critical number, the three-body interaction takes control 
and the number of atoms suddenly drops below the crit- 
ical level by recombination loss signaling a collapse 
Then the absorptive term takes over and the number of 
atoms starts to increase again. This continues indefi- 
nitely showing an infinite sequence of collapse. 



II. COUPLED GROSS-PITAEVSKII EQUATION 
WITH ABSORPTION 

We consider the following spherically symmetric cou- 
pled GP equation with two components at time t for the 
condensate wave function (r, r) 



h 1 d'^ 1 2 2 TV r I / / M 2 

-ih^ 4>,{r,T)=0, (2.1) 



i — 1,2, where m is the atomic mass. Here gij = 
Airli^ Gij / m is the coupling constant for atomic interac- 
tion, Nj the number of condensed atoms in state j, and 
oj the frequency of the harmonic oscillator trap. The pa- 
rameter Ci has been introduced to modify the frequency 
of the trap for the atoms in each quantum state. 

As in Refs. g it is convenient to use dimcnsionless 
variables defined by x = V^t/qy^q , and t — tlo, where 
ajjQ = ^JyTjJjruJ), and (^i{x,t) = x?/'i(?-, r)(V27ra^^)i/2 . 
In terms of these variables Eq. (|2.1|) becomes pi 
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\^j{x,t)f^ 



|0,(x,i)|^ 



dt 



(j)i{x,t) = 0, 



(2.2) 



where = 2^/2N jUij / aY^^ could be negative (positive) 
when the corr esp onding interaction is attractive (repul- 
sive) . In Eq. (^j) we have introduced a diagonal absorp- 
tive i7i and a quartic three-body term —\£^i\(t)i{x^t)\'^/x'^ 
appropriate to study collapse g]. For 7i = = 0,i = 
1,2, the normalization condition of the wave function is 



|0,(x,t)|2dx = 1. 



(2.3) 



The root-mean-square (rms) radius of the component i 
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(2.4) 



III. NUMERICAL RESULTS 



To solve Eq. (2.2) we discretize it in both space (us- 
ing step 0.0001) and time (using step 0.05) employing a 
Crank-Nicholson-type rule and reduce it to a set of al- 
gebraic equations which is then solved by iteration using 
the known boundary conditions, e.g., 10^(0, i)| — 0, and 
Vaxix^oo |0i(a^,i)| ^ exp(— a;^/4). The iteration is started 
with the known normalized (harmonic oscillator) solu- 
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{t) at time t is defined by 



tion of Eq. ( |2.2| ) obtained with = at t = 0. The 
nonlinear constants in this equation are increased by 
equal amounts over 500 to 1000 time iterations starting 
from zero until the desired final values are reached. This 
iterative method is similar to one in the uncoupled case 
[^,|| . A detailed account of the numerical procedure for 
the coupled case will appear elsewhere. 



A. Stationary Problem 

First we consider the stationary solution of Eq. ( |2.2| ) 
with 7i = = 0, which illustrate the collapse. As the 
three scattering lengths and two numbers Ni are all 
independent, the four parameters riij are also so with one 
restriction: the signs of ni2 and n2i are identical. 

Now we study the simplest case of collapse by taking 
only the interaction between the atoms in state 1 to be 
attractive corresponding to a negative an . All other scat- 
tering lengths — 022 and ai2 (= 021) — are taken to be 
positive. Quite expectedly, here the first component of 
the condensate could experience collapse. Although the 
present formulation is generally valid, one has to choose 
numerical values of the parameters before an actual cal- 
culation. 

The collapse of the first component is illustrated in Fig. 
1 (a) for mi = —3.814, 7122 = 4 ni2 = "21 = 1, ci = 0.25, 
C2 = 4. These parameters are in dimcnsionless units and 
one can associate them with an actual physical problem 
of experimental interest. For this we consider the state 1 
to be the states of ''Li with attractive interaction as in the 
actual collapse experiment with |aii|/ajj^Q — 0.0005 
As nil ~ 2-\/2A^i|aii|/aj^Q this corresponds to a boson 
number iVi ~ 2700. This number is larger than the max- 
imum number atoms permitted in the BEG of a single 
component ^Li which is about 1400 ||^. The presence of 
the second component with repulsive interaction allows 
for a formation of a stable BEG with more ^Li atoms in 
quantum state 1 than allowed in the single-component 
BEG. Similar conclusion was reached by Esry |]l3| in a 
study of a coupled BEG in a different context. We find 
from Fig. 1 (a) that (pi is very much centrally peaked 
compared to (j)2- This corresponds to a small rms radius 
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and large central density for (pi denoting an approxima- 
tion to collapse. If the number Ni is slightly increased 
beyond 2700 the first component of the condensate wave 
function becomes singular at the origin and no stable 
stationary solution to Eq. (2.2) could be obtained. 

Next we discuss the collapse by taking only the interac- 
tion among atoms in two different states to be attractive 
corresponding to a negative 012 (= 021). The atomic 
interaction in both quantum states 1 and 2 is taken to 
be repulsive corresponding to a positive an and 022. Al- 
though it is a problem of theoretical interest for the study 
of collapse, it has no experimental analogue in terms of 
^Li. We illustrate the approximation to collapse in this 
case in Fig. 1 (b) for parameters nu = 1,7122 = 1-5, 
ni2 = —5.95,7121 = —2, ci = 1, C2 = 0.25. Both 
wave-function components are peaked near a; = and 
have small rms radii. The system would collapse with 
a small increase of |ni2| and/or |?t.2i|. Here the interac- 
tions among atoms in states 1 and 2 are both repulsive. 
The collapse is a consequence of the attraction between 
an atom in state 1 and one in state 2. This leads to 
a dominance of nonlinear off-diagonal coupling terms in 
the coupled GP equation. 

Finally, in Fig. 1 (c) we illustrate the approximation to 
collapse of both components when all scattering lengths 
are negative. This corresponds to taking all possible in- 
teractions attractive. The parameters in this case are 
nil = = -1, 7112 = "21 = -0.552, ci = 4, C2 = 0.25. 
This has an experimental analogue in terms of two states 
of ^Li. We assume the atomic interaction in both states 
to be equally attractive corresponding to a negative scat- 
tering length: an = a22- For |aii|/ajjQ ~ 0.0005 as in 
the actual experiment Q, one has Ni = N2 — 700. The 
total number of particles in this case is roughly 1400, 
which is equal to the critical number observed in the ac- 
tual experiment of collapse in ''Li. Both wave-function 
components could become singular in this case as all pos- 
sible interactions are attractive. 



B. Time-dependent Problem 

Although the collapse of the coupled condensates could 
be inferred from the shape of the stationary wave func- 
tions of Fig. 1 (sharply peaked centrally with small rms 
radii), we also study the dynamics of collapse from a time 
evolution of the full GP equation (2.2) in the presence 
of an absorption and three-body recombination, e.g., for 
7i ^ and 7^ as in the uncoupled case For this 
purpose we consider the solution of Eq. (2^) normahzed 
according to Eq. (2.3) at t ~ obtained with ji ~ £,i = 
and allow this solution to evolve in t ime with 7^ 7^ and 

7^ by iterating the GP equation (2.2). The fractional 
change in the number of atoms due to the combined ef- 
fect of absorption and three-body recombination is given 
by 
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and the rms radii by Eq. (2.4). The continued growth 
and decay of the number of particles in the condensate 
would signal the possible collapse in a particular case. 
The oscillation of the rms radius would demonstrate the 
consequent radial vibration of the condensate. 

Now we study the time evolution of the number of 
atoms of the two components and the corresponding rms 
radii. The general nature of time evolution is indepen- 
dent of the actual values of 7^ and employed provided 
that a very small value for 0.001) and a relatively 
larger one for jil^ 0.01 to 0.1) are chosen The follow- 
ing parameters were chosen in case of models (a), (b), and 
(c) of Fig. 1: (a) and (b) 71 = 72 = 0.03,6 = 6 = 0.001, 
(c) 71 = 0.15,72 = 0.03,6 = 0.002,6 = 0.003. The 
fractional change in the number of atoms for the two 
components are shown in Figs. 2 (a), (b), and (c). The 
results for <t < 100 in Fig. 2 are calculated with 2000 
iterations of the GP equation (2.2) using a time step 0.05. 

The quadratic nonlinear terms in model (a) are all re- 
pulsive in channel 2, the corresponding wave function 
(02 ) of Fig. 1 (a) does not show any sign of approxi- 
mation to collapse as in channel 1 where the diagonal 
nonlinear term is attractive. The results reported in Fig. 
2 (a) are consistent with this. The number of particles 
A^i of the first component undergoes successive growth 
and decay, whereas that of the second component keeps 
on growing indefinitely typical to a repulsive interaction. 

For model (b) the effective nonlinear terms in channels 
1 and 2 are both repulsive and it should be possible to 
have collapse in both channels by decreasing the dom- 
inating off-diagonal quadratic nonlinear terms 7112 and 
7121 corresponding to an increase in attraction between 
an atom in state 1 and one in state 2. However, for the 
actual parameters of this model only the component 1 
exhibits collapse. This is consistent with the more sin- 
gular nature of 4>i reported in Fig. 1 (b), compared to 
(j)2- Consequently, in Fig. 2(b) only component 1 ex- 
periences collapse; the number of particles N2 keeps on 
growing with time. 

In model (c) all the quadratic nonlinear terms are at- 
tractive. Consequently, in Fig. 2 (c) we find a series of 
collapse in both channels. The collapse is most favored 
in model (c) with attractive diagonal and nondiagonal 
nonlinear terms. This corresponds to attraction between 
two atoms in state 1, between two atoms in state 2, and 
between an atom in state 1 and another in state 2. The 
next favored case is of model (a) where the diagonal non- 
linear term is negative in channel 1. Here only the atomic 
interaction in state 1 is attractive, all other atomic inter- 
actions are repulsive. The least favored case is of model 
(b) where only the off-diagonal nonlinear terms are nega- 
tive. This corresponds to repulsion between two atoms in 
state 1, and between two atoms in state 2, and attraction 
between an atom in state 1 and another in state 2. In 
the last case, collapse takes place due to the dominance 
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of the attractive nondiagonal nonlinear term over the re- 
pulsive diagonal one in channel 1. This is explicit in Fig. 
2 where the frequency of collapse decreases from model 
(c) to (a) and then to (b). 

Finally, in Figs. 3 (a), (b), and (c) the rms radii for 
the two components are shown for models of Figs. 2 (a), 
(b), and (c), respectively. In case of models (a) and (b) 
we find from Figs. 2 (a) and (b) that the number N2 
grows with time. This is reflected in the growth of the 
corresponding rms radii in Figs. 3 (a) and (b). In case of 
model (c) there is collapse in both channels and both the 
rms radii oscillate with time. This radial vibration of the 
collapsing condensate(s) also takes place in the uncoupled 
case However, from Figs. 3 (a) and (b) we find that 
due to a collapse in one of the channels, both rms radii 
could execute oscillations. In one of the channels it is a 
direct consequence of collapse, in the other it is due to a 
coupling to the channel experiencing collapse. 

IV. CONCLUSION 

To conclude, we studied the collapse in a trapped BEC 
of atoms in states 1 and 2 using the GP equation when 
some of the atomic interactions are attractive. We mo- 
tivate parts of this study with two atomic states of ^Li. 
The component i of the condensate could experience col- 
lapse when the interaction among atoms in state i is 
attractive. Both components could experience collapse 
when at least the interaction between an atom in state 1 
and one in state 2 is attractive. The collapse is predicted 
from a stationary solution of the GP equation. The time 
evolution of collapse is studied via the time-dependent 
GP equation with absorption and three-body recombi- 
nation. The number of particles of the component (s) of 
BEC experiencing collapse alternately grows and decays 
with time. With the possibility of observation of coupled 
BEC, the results of this study could be verified experi- 
mentally in the future. 
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Figure Caption: 

1. Wave function components 4>i{x) (full line) and 
4>2{x) (dashed line) vs. x for two coupled GP equations 
with (a) nil = —3.814, 7122 — 4, ni2 = ^21 = 1, ci = 0.25, 
C2 = 4; (b) nil = 1,^22 = 1-5, ni2 = -5.95, n2i = -2, 
Ci = 1, C2 = 0.25; and (c) nn = 7222 = -1, ni2 = n2i = 
-0.552, ci = 4, C2 = 0.25. 

2. The fractional change in the number of atoms 
Ni{t)/Ni{0) vs. t for component 1 (full line) and 2 
(dashed line) for models (a) and (b) with 71 = 72 = 0.03 
and Ci = ^2 = 0.001, and for (c) with 71 = 0.15,72 = 
0.03, ^1 = 0.002, and ^ = 0.003. The parameters are as 
in Fig. 1. 

3. The time dependence of rms radii XY[n.g{t) of models 
(a) . (b) . and (c) for component 1 (full line) and 2 (dashed 
line). The parameters are as in Figs. 1 and 2. 
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